Search for the Rare Decay Kl — > 7r°7r°7 
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The KTeV E799 experiment has conducted a search for the rare decay Kl — » 7r°7r°7 via the 
topology Kl — > tv°-k%'j (where tt% —* r ye + e~). Due to Bose statistics of the 7r° pair and the real 
nature of the photon, the Kl — * •K°ir°y decay is restricted to proceed at lowest order by the CP 
conserving direct emission (DE) of an E2 electric quadrupole photon. The rate of this decay is 
interesting theoretically since chiral perturbation theory predicts that this process vanishes at level 
0(p 4 ). Therefore, this mode probes chiral perturbation theory at 0(p 6 ). In this paper we report a 
determination of an upper limit of 2.43 x 10~ 7 (90% CL) for Kl — * Tv°-K 0r y. This is approximately 
a factor of 20 lower than previous results. 

PACS numbers: 13.20.Eb, 13.25.Es 



I. INTRODUCTION 



While measurements of branching ratios for kaon deca' 
modes such as Kl — > 77 have shown good agreement [1| 
with chiral perturbation theory (ChPT) calculations to 
order 0(p 4 ), ChPT calculation to order 0(p 6 ) have been 
difficult to check due to the presence of the large lower 
order terms. Since the Kl — ► 7r°7r°7 mode is forbidden 
to order 0(p 4 ) in ChPT 0, this mode presents an op- 
portunity to check ChPT at higher orders. Therefore 
an experimental measurement of this mode's branching 
ratio is of interest. At present, a chiral dimensional anal- 
ysis Q for the Kl tt°tt°j mode results in a branching 
ratio of 7 x 10~ n . Because of the Bose statistics of the 
7r°7r° pair and the presence of a real direct emission 7, the 
Kl — > 7r°7r°7 decay proceeds at lowest order via the CP 
conserving direct emission of a electric quadrupole (E2) 
photon. This causes a large suppression of this mode rel- 
ative to the associated Kl — > 7r + 7r~7 mode. Accordingly, 
observation of Kl — > 7r°7r°7 would give a measurement of 
E2 quadrupole emission which is difficult to extract from 
the El and Ml dominated Kl — » tt + it~j mode. The 
Ml amplitude determined from Kl —* Tr + n~j decay to- 



gether with a plausible estimate for the E2 amplitude 
given by 

gE2 (pi -P2) ■ k 
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can be used to obtain an estimate of the branching ratio 
for Kl — > 7r 7r°7. Here p\ and P2 are the momenta of the 
pions and k and e are the momentum and polarization 
of the photon respectively, qei is the coupling of the E2 
amplitude. A -1 , a mass parameter that measures the ex- 
tent of the region over which the interaction takes place, 
is assumed to be of order of the p meson mass. The CP 
violating direct emission of a M2 magnetic quadrupole 
photon is much smaller than the E2 transition and can 
be ignored. According to Ref. the decay rate for 
Kl — * 7r°7r°7 can be estimated by comparing it to the 
Ml direct emission (DE) amplitude of the Kl — > 
decay 



(2) 



Employing this amplitude, the ratio 
T(K L -> 7r°7r°7) 
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can be formed. Using the fraction of T(Kl — » vr + 7r 7) 
to Y(Kl — > 7T + 7r _ ) due to Ml direct emission of (14.2 ± 
0.28) x 1(T 3 from Ref. % and a K L -> tt+tt- branch- 
ing ratio of (2.090 ± 0.025) x 10~ 3 from the PDG @, a 
branching ratio for Ml DE part of the Kl — > 7r + 7r~7 
decay of 2.96 x 10~ 5 {E 1 > 20MeV) is obtained. Using 
Eq. 3 above, this gives an estimated branching ratio of 
1.08 x 10 -8 for Kl — > 7r°7r 7 (under the assumption that 
gE2 has similar magnitude to gun)- 

The KTeV collaboration previously reported an upper 
limit of of 5.4 x 10~ 9 (90% CL) 7] for the related decay 
Kl — » 7r°7r°e + e~ in which the direct emission E2 photon 
was virtual. However, we point out that this decay, in 
contrast to the Kl — > 7r°7r°7 decay where the photon is 
real, can proceed via an additional process forbidden for 
a real photon. The virtual photon is emitted in a J=0 
state from the Kl allowing a transition of the Kl to a 
Ks followed by the CP allowed decay K§ — > 7r°7r°. This 
is the so-called charge radius amplitude 0]. The rates 
for these two decays would have similar branching ratios 
except for the the charge radius amplitude. 

The previous best upper limit of 5 x 10~ 6 (90% CL) 
for the Kl — > 7r°7r°7 decay obtained by the NA31 ex- 
periment 0] was achieved by searching for this mode in 
events with five photons in the final state. Due to the 
particular configuration of triggers and prescale factors 
in the KTeV experiment, the most sensitive method was 
to require that one of the neutral pions undergo Dalitz 
decay ir° — > e + e~~f. 

II. THE KTEV E799 EXPERIMENT 

The search for the Kl — > Tv Tr L \'y mode was performed 
using the 1997 and 1999 runs of KTeV E799 II at Fermi 
National Accelerator Laboratory. The KTeV E799 ex- 
periment used two almost parallel Kl beam lines gener- 
ated by interactions of 800 GeV/c protons in a BeO tar- 
get some 90 meters upstream of the KTeV spectrometer. 
The KTeV spectrometer consisted of a 70 meter evacu- 
ated decay tube followed by, in sequence, two stations 
of drift chambers, a large aperture dipole magnet, two 
more stations of drift chambers, a multi-plane transition 
radiation detector, a 3100 element Csl electromagnetic 
calorimeter and a muon detector. 

The 1997 and 1999 runs differed in the following ways. 
The spill length was doubled from 20 in the 1997 run to 
40 seconds in the 1999 run. The proton intensity on the 
BeO neutral kaon production target was increased from 
4 x 10 12 in the 1997 run to 6 - 10 x 10 12 per spill in the 
1999 run. Another important difference between the 1997 
and 1999 run was that the magnetic field was decreased 
from a transverse deflection of 205 MeV/c in 1997 to 150 
MeV/c in 1999 to increase acceptance for some of the 
neutral kaon modes. The trigger which demanded two 
charged tracks and at least four electromagnetic clusters 



was loosened midway during the 1997 run by changing 
the minimum thresholds for energy in the electromag- 
netic calorimeter and remained loose for the 1999 run. 

II A. The K L -» tt°tt%j Signal Criteria 

The Kl — > 7r°7r^,7 final state consists of four pho- 
tons plus an e + e~ pair. The KTeV analysis required 
six electromagnetic showers (E > 0.6 GeV) in the Csl 
calorimeter, two of which were associated with two re- 
constructed charged tracks. The two charged tracks were 
required to be consistent with the two track trigger, form 
a good charged track vertex, have opposite charge, and 
have 0.95 < E/p < 1.05 (where E is the energy of the 
associated shower in the Csl calorimeter and p is the mo- 
mentum determined from magnetic bending). 

Events passing these cuts were required to satisfy sev- 
eral additional cuts to select Kl — * 7r°7r°7 and to re- 
ject the major background due to Kl — + 7r 7r 7r2, de- 
cays. These events can contribute to the background if 
one of the photons was not detected or two of the pho- 
tons reconstructed as single photon. Additional cuts were 
performed including a vertex cut determined by cycling 
over all 77 and e + e~7 combinations in a given event 
for all possible vertex positions to determine the vertex 
where the M e + e - 7 and M 77 were closest to the ir° mass. 
The vertex thus obtained was required to be between 95 
and 150 meters from the BeO target in the KTeV de- 
cay volume and to have a good vertex \ 2 ■ Another cut 
was made on the variable {Pl)-k° defined as the longi- 
tudinal momentum squared of the 7r° in the frame in 
which the momentum of the 7r 7r^ pair is totally trans- 
verse. If the event is a Kl — * 7r 7r 7r^ decay, {Pl)-k° IS 
greater than zero. If the event is a Kl — > 7]"°7]£>7 decay, 
(P 2 ) w o is less than zero (modulo resolution). Accordingly 
(P|) w o was required < —0.005 GeV 2 /c 2 in order to select 
Kl — » 7r°7r^ ) 7 decays and reject Kl — * 7r 7r 7r^,. 

The analysis cuts were varied slightly between the 
1997 and 1999 run. The main changes were to alter the 
shape of the p 2 (eejjjj) and M ee7777 signal region and 
to tighten the decay vertex position and the overlapping 
shower cuts to reject more background. 

II B. Backgrounds 

Monte Carlos were developed for the signal Kl —> 
7r ° 7r r)7 m ode incorporating the model of Ref. Q as well 
as for the various backgrounds, and the normalization 
mode incorporating all known features of the E799 spec- 
trometer and beam. The various Kl modes were gener- 
ated using experimentally determined decay parameters. 
Trigger and analysis cuts applied to the Monte Carlos 
were the same as those applied to the data. Background 
and normalization mode Monte Carlos were generated to 
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FIG. 1: a) comparison of cos(#) for a Kl — * 7r 7r^,7 signal 
MC vs. cos(6 l ) as observed in the 1997 data; b) comparison 
of cos(#) for the Kl — > 7r°7r°7ro background MC vs cos(#) for 
the 1997 data 



gated to see if there was any contribution to background 
after the cuts to select K L — > 7r 7r^,7. These modes 
included K L — » ir^n^n^ where a Dalitz pair is unde- 
tected, if £ — * 7r 7r 7r° followed by the decay of a 7r° into 
an e + e~ pair, i^L — * 7r°7r 7r° followed by the decay of 
a 7T° into e + e~e + e~ with one of the e + e~ pairs unde- 
tected, K L — > 7r°e + e _ 7 decays with accidental photons, 
and ifi — > 7r°7r^ plus accidental photons. In no case did 
any of these modes contribute significantly to the back- 
ground. The agreement of the Monte Carlo simulation of 
the remaining background with the 1997 and 1999 data 
outside the signal region is shown in Fig. [2] 



III. SEARCH FOR A K L -> Tr ^ SIGNAL 



The signal regions for the 1997 and 1999 data were 



based on the M, 



ee7777 and pi (ee7777) resolutions calcu- 
lated using a signal mode Monte Carlo. Here pf is mea- 
sured relative to the direction of the Kl determined by 
the line connecting the BeO target center and the decay 
vertex. For the 1997 data, the signal region was chosen to 
be rectangular with 0.494 GeV/c 2 < M < 0.501 GeV/c 2 
and p\ < 0.00015 (GeV/c) 2 . In the later analysis of the 
1999 data, it was decided to use a contour containing 68% 
of the signal Monte Carlo events as determined from a 
joint probability distribution based on the signal Monte 
Carlo resolutions for p^(eejjjj) and Me e7 -v 77 . Since the 



ee7777 . 

signal box had already been opened for the 1997 data 
(and no events were observed as discussed below), the 
rectangular signal region was kept for the 97 data. 

As shown in Fig. [3^, no events were observed in the 
1997 data in this region when the box was opened af- 
ter the blind analysis was completed. From the K L —> 
7r 7r 7r£, background Monte Carlo, we expected 0.83±0.41 
background events in the 1997 run. As shown in Fig.[3jD, 
one data event was found in the 1999 signal region while 
the background K L — ► 7r 7r°7r2j Monte Carlo yielded no 
events in the signal region for four equivalent K L fluxes. 
These Monte Carlo background estimates for the 1997 
and 1999 runs are consistent with the background esti- 
mates based on mass and p\ sideband projection from 
the data. 



■+7r°7r° 7 
7r u 7r u 7r^ with the 



produce Monte Carlo samples several times the 1997 and 
1999 run K L fluxes. 

As indicated in Fig.[TJ a comparisons of a K L 
and a background Monte Carlo K L 
1997 data shows that the angle 8 between the direct emis- 
sion photon and one of the two 7r°'s in the 7r 7r° center of 
mass is peaked for the signal mode at cos# = ±0.6 and 
flat for K L -> 7r°7r°7r^. The same is true for the 1999 
data. 

Therefore a final cut to select Kl — * 7r°7T£,7 was made 
requiring events to have 0.25 < |cos6*| < 0.9 to improve 
signal to background. Several other decays were investi- 



IV. DETERMINATION OF THE UPPER LIMIT 

FOR K L 



7^7 



Since the Kl — ► 7r 7r 7r^, branching ratio is well mea- 
sured and since events where a photon was lost down 
the beam hole had the same topology as the signal and 
were copious, these events were used for the normaliza- 
tion for the experiment. The reason for selecting this 
topology was to minimize the systematics due to trigger 
and event configuration between the signal and normal- 
ization modes. 
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FIG. 2: Data and Monte Carlo background comparisons pass- 
ing all cuts except the signal region contours of Fig. [3] for a) 
the 1997 run; b) the 1999 run. 



FIG. 3: 

b) M ee ~ 



a) M e , 



vs. pt (ee7777) for the 1997 KTeV data; 



vs. p t (ee7777) for the 1999 KTeV data 



The criteria to select these normalization mode events 
included the same topology criteria, the vertex position 
cut, and the tt° mass cuts as the signal mode. Cuts that 
were special to the selection of the normalization mode 
included a requirement that (P|) ff o > 0.0. In addition, 
to insure that the event was consistent with a Kl —> 
Tr ^ ^ decay with a photon going down the beam hole, 
the missing photon was reconstructed by assuming a kaon 
mass and net zero pi for the event. This resulted in two 
solutions for energy and trajectory of the missing photon. 
If one of the solutions for the missing photon combined 
with the unpaired photon in the observed event to form 
a 7T° mass within 0.05 GeV/c 2 , if the other three photons 
and the e + e~ in the event paired to form 7r°'s, and if the 
missing photon's calculated trajectory pointed to cither 
beam hole, the event was included in the normalization 
sample. Normalization mode yields of 7725 and 9024 
K L -> 7r°7r°7r£, were obtained for the 1997 and 1999 data 
using these criteria. 
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IV A. The Single Event Sensitivity for K L 
Signal 



We can avoid the extra error due to the branching 
ratio uncertainty in the Dalitz decay by noting that the 
Kl — > 7T°7r^7 single event sensitivities (SES) are given 
by 



SES(n\° Dl ) 97 , g9 



1 



A(n0n° Dl ) 97t99 N(K 



Lj97,99 



(4) 



where A(7r 7r^,7)97 i 9g are the acceptances for the Kl — > 
7r ° 7r £)7 Dalitz decay mode and N(-Ki)97,9g are the Kl 
fluxes for the 1997 and 1999 modes. The fluxes are 



N(3k° d ) 



norm 
97,99 



(5) 



so the SES's for K L 



7r°7r°7 



SES. 



97,99 



WR(K L 



77777 are given by 
\ tt )BR{ir a -> 77)A(3tt 



.CM97.99 



2A( 7 r°7r0 ) 7) 97 , 99 iV(37r0 3 ) 9 l 7 o T m 



',99 



(6) 

where N(3n%)$fi$ and A(37r^) 9 l f' 9 ' 9 t are the number and 
the acceptance of normalization mode iff, — > 7r 7r 7r£, 
events in the 97 and 99 runs. The factor of 2 and 3 arise 
from the number of 7r°'s that can produce a Dalitz pair 
in the signal and normalization modes. Note that the 
Dalitz decay branching ratio cancels out in the SES for 
K L -> 7r°7r°7. 



From the Kl 



7^2,7 



and Kl 



' D 



Monte 



Carlos, the signal and normalization acceptances were 
determined to be 0.101% and 4.21 x 10~ 6 for the 1997 
run and 0.085% and 3.26 x 10~ 6 for the 1999 runs re- 
spectively. Using Eq. 6 and the measured 19.56±0.12% 
branching ratio for K L -> 7r 7r°7r° 10], the 1997 and 1999 
SES's were determined to be 1.56 x 10~ 7 and 1.25 x 10~ 7 
respectively. In addition, Kl fluxes of 2.68 xlO 11 and 
3.99 x 10 11 were obtained for the 1997 and 1999 runs. The 
SES's for 97/99 were combined using 
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1 



Error on the 3tt branching ratio 
Error on the ir° — » 77 branching ratio 
Normalization Monte Carlo/data disagreements 
Signal/normalization Monte Carlo disagreements 



Total Systematic Error 



0.61% 
3.24% 
3.57% 
5.35% 



7.23% 



TABLE I: Systematic uncertainties for K L — > 7r°7r°7 SES 



Monte Carlo. A disagreement in signal mode data and 
Monte Carlo could not be checked directly since no signal 
is observed. However, from an inspection of the calcula- 
tion of the SES, it can be seen that the SES is propor- 
tional to (N signal /A signal )/(N norm /A norm ). So, to the 
level that the normalization and signal mode topologies 
are similar, it was expected that differences in acceptance 
between the Monte Carlo and the data would tend to 
cancel. Studies were done of the effect of disagreements 
between the signal mode and normalization mode Monte 
Carlos by adjusting the Monte Carlos until they agreed 
and determining the effect on the result. 

The percentage systematic errors in the SES due to 
these sources are shown in Table |U Adding these sys- 
tematic errors in Table I in quadrature, we obtain a total 
systematic error in the SES of 7.23%. Taking into ac- 
count the statistical error of 1.09%, this leads to a total 
error of 7.31% resulting in a SES error of 5.07 x 10~ 9 . 



IV C. Calculation of the Upper Limit for 

Kl -> 7r°7r£,7 Signal 



The method of Ref. 11| was used to obtain an upper 
limit for Kl — > 7r°7r°7. If n exp is the expected number of 
signal plus background events in the signal box and ribkg 
is the expected number of background events, then the 
probability for observing n events is given by a Poisson 



SES, 



combined 



SES 



97 



SESq 



(7) distribution P(n exp ,n) whose mean is 



resulting in a combined SES of 6.93 x 10 



nbkg 



BR(K L -► 7r 7r°7)) /SES (8) 



IV B. Systematic Errors 

Systematic errors for the upper limit could arise from 
several sources. Among them was the error in the branch- 
ing ratio for Kl —* 7r 7r°7r°. There were also disagree- 
ments between the distributions of the Kl — » ir ^^ 
normalization mode data and Monte Carlo. The effect of 
these disagreements on the SES were studied by adjusting 
the Monte Carlos to eliminate the disagreements and see- 
ing what changes in the SES took place. There could also 
be disagreements between the Kl — > ^ q ^l>1 data and its 



The expected and observed backgrounds for the 97 and 
99 runs were added to form the final observed and 
expected backgrounds. A confidence region was con- 
structed containing 90% of the Poisson distribution by 
varying the Kl — > 7r°7r°7 branching ratio between to 
5 x 10~ 7 . Using this confidence region, a 90% CL upper 
limit for the BR(K L -> 7r°7r°7) of 2.43 x 10~ 7 was de- 
termined, assuming the E2 model of Ref. The SES 
error was incorporated by varying the SES over the 90% 
CL range of the SES indicated by its combined statistical 
and systematic error. 
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V. CONCLUSIONS 

In conclusion, the KTeV collaboration has obtained 
an upper limit for the branching ratio for the rare decay 
K L -> 7r°7r^7 of 2.43 x 1(T 7 using the 1997 and 1999 data 
and assuming that the decay proceeds mainly via direct 
E2 photon emission. This limit is approximately twenty 
times lower than the best published upper limit Q, pro- 
viding a much more stringent upper limit for ChPT the- 
oretical calculations of the E2 amplitude. 
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